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1 Introduction 

Dark matter (DM) is estimated to account for about 23% of the total mass of the universe, and 
to be five times more abundant than the known baryonic matter. While the existence of DM 
is inferred from astrophysical observations, there is very little information about its nature or 
how it interacts with ordinary matter. 

In this paper, we consider a simplified scenario QJjSl in which DM has a particle explanation 
and, in particular, there is only one new Dirac fermion related to DM within the energy reach 
of the LHC. The fermion interacts with quarks via a four-fermion contact interaction, which 
can be described by an effective field theory (EFT) Lagrangian: 

^ = EEC,« Wq)(*Tfr), (1) 

<7 i 

where C represents the coupling constant, which usually depends on the scale of the interac¬ 
tion (M*). The operator T describes the type of the interaction, including scalar (T = 1), pseu¬ 
doscalar (T = y 5 ), vector (T = y"), axial vector (T = y^y 5 ), and tensor interactions (T = a^ v ). 
The exact value of the constant C depends on the particular type of the interaction. 

This scenario can lead to the production of DM particles in association with a hard parton, a 
photon, or a W or Z boson. The first two production modes are usually referred to as mono¬ 
jets |Tjj3«6] and monophotons JU, respectively. Recent monojet results from the ATLAS (Zl and 
CMS HI Collaborations have placed lower limits on M* for some typical couplings in Eq. ([lj. 
The ATLAS Collaboration @ has also searched for DM particles in events with a hadroni- 
cally decaying W or Z boson. Assuming a DM particle with a mass of 100 GeV, the excluded 
interaction scales are below about 60 GeV |9|, 1040 GeV (U, 1010 GeV 03 , and 2400 [j9| GeV for 
scalar, vector, axial-vector, and tensor interactions, respectively, and the excluded scale is below 
410 GeV [[8] for a scalar interaction between DM particles and gluons. 

The exclusion limit for a scalar interaction between DM particles and quarks is the least strin¬ 
gent among all the interaction types that have been probed. In this interaction the coupling 
strength is proportional to the mass of the quark: 

7YI 

L mt = (2) 

As a consequence, couplings to light quarks are suppressed. A recent paper IflOl suggested 
that the sensitivity to the scalar interaction can be improved by searching in final states with 
third-generation quarks. It has also been noted that the inclusion of heavy quark loops in the 
calculation of monojet production IfTTI increases the expected sensitivity. 

In this paper, we report on a search for the production of DM particles in association with a pair 
of top quarks, and consider only the scalar interaction. The ATLAS Collaboration has recently 
searched for DM particles in association with heavy quarks fl2l , placing more stringent limits 
on the scalar interaction between DM particles and quarks than the mono-W/Z search |9|. 
Assuming a DM particle with a mass of 100 GeV, the excluded interaction scale is 120 GeV for 
scalar interaction between top quarks and DM particles. Figure [T] shows the dominant diagram 
for this production at the LHC. In this paper we focus our search on events with one lepton 
(electron or muon) in the final state. 
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3 Data and simulated samples 



Figure 1: Dominant diagram contributing to the production of DM particles in association with 
top quarks at the LHC. 

2 The CMS Detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintilla¬ 
tor hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Extensive forward calorimetry complements the coverage provided by the barrel 
and endcap detectors. A more detailed description of the CMS detector, together with a def¬ 
inition of the coordinate system used and the relevant kinematic variables, can be found in 
Ref. HU. 

3 Data and simulated samples 

The data used in this search were recorded with the CMS detector at the LHC at i/s = 8 TeV, 
and correspond to an integrated luminosity of 19.7 fb 1 . The data were collected using single¬ 
electron and single-muon triggers, with transverse momentum (pj) thresholds of 27 and 24 GeV, 
respectively. The efficiencies of these triggers in data and simulation are compared, measured 
using a tag-and-probe method [T4l . and correction factors are applied to the simulation. 

DM signals are generated with MADGRAPH v5.1.5.11 Ifl5l leading order (LO) matrix element 
generator using the CTEQ6L1 parton distribution functions (PDF) fl6l . The dominant standard 
model (SM) background processes for this search are tt+jets, tt + y/W/Z, W+jets, single top 
quark, diboson (WW, WZ, and ZZ) and Drell-Yan events. All of these backgrounds except 
single top quark and WW events, are generated with the MADGRAPH using CTEQ6L1 PDF. 
The top-quark pj distributions in the tt+jet sample generated from MADGRAPH are reweighted 
to match the CMS measurements, following the method described in Ref. IH7I . Single top quark 
processes are generated with the next-to-LO (NLO) generator POWHEG vl.O using the CTEQ6M 
PDF HT6| . The WW background is generated with the PYTHIA v6.424 I flBl . All events generated 
with MADGRAPH are matched to the PYTHIA lfl8| parton shower description. All events are 
passed through the detailed simulation of the CMS detector based on GEANT4 v9.4 Iil9l . 

The cross sections of tt+jets Il20l and W/Z+jets ||2H backgrounds are calculated at next-to-NLO. 
Other backgrounds are calculated at NLO. The single top quark cross section is taken from 
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Ref. It22l . the tt + Z cross section from Ref. Il23l . the ft + W cross section from Ref. [24ll , the tt + 7 
cross section from Ref. lf25H and the diboson cross sections are from Ref. l26l . 

Additional minimum bias events in the same LHC bunch crossing (pileup) are added to all 
simulated events, with a distribution in number matching that observed in data. 


4 Object reconstruction 

A particle-flow (PF) based event reconstruction Ii27l l28fl is used by CMS, which takes into ac¬ 
count information from all subdetectors, including charged-particle tracks from the tracking 
system and deposited energy from the ECAL and HCAL. Given this information, all parti¬ 
cles in the event are classified into mutually exclusive categories: electrons, muons, photons, 
charged hadrons, and neutral hadrons. Primary vertices are reconstructed using a determinis¬ 
tic annealing filter algorithm [29), with the event primary vertex defined as the vertex with the 
largest sum of the squares of the pj of the tracks associated with that vertex. 

Electron candidates are reconstructed from energy clusters in the ECAL matched with tracks 
Mi- The electron trajectory in the tracker volume is reconstructed with a Gaussian sum fil¬ 
ter [3D algorithm that takes into account the possible emission of bremsstrahlung photons 
in the silicon tracker. The electron momentum is then determined from the combination of 
ECAL and tracker measurements. Electrons are identified by placing requirements on the 
ECAL shower shape, the matching between the tracker and the ECAL, the relative energy 
fraction deposited in HCAL and ECAL, the transverse and longitudinal impact parameters 
of the tracker track with respect to the event primary vertex, photon conversion rejection, and 
the isolation variable Rf so . The isolation variable is defined as the ratio to the electron trans¬ 
verse momentum, of the sum of pj of all other PF candidates reconstructed in a cone of radius 
A R = \/(A/ 7) 2 + (A cp) 2 = 0.3 around the electron candidate, where // is the pseudorapidity 
and cf) is the azimuthal angle. The pj sum in the isolation cone is corrected for the contribu¬ 
tions of pileup interactions on an event-by-event basis. Isolated electrons satisfy Rf so < 0.1. 
The electron is required not to be in the transition region between the barrel and the endcap 
ECAL (1.44 < \rj\ < 1.57) because the reconstruction of an electron object in this region is not 
optimal [30) . After all these requirements, electrons are selected if they satisfy pj > 30 GeV and 
\v\ < 2-5. 

Muon candidates are reconstructed by combining tracks from the tracker and muon system [32), 
resulting in "global-muon tracks". The PF muons are selected among reconstructed muon 
track candidates by imposing minimal requirements on the track components in the muon sys¬ 
tem and taking into account matching with small energy deposits in the calorimeters Il27l l28l. 
Muons from cosmic rays and from light hadrons that decay in flight, or from b hadrons, and 
hadrons misidentified as muons are suppressed by applying requirements on the quality of the 
global-muon fit, the number of hits in the muon detector and in the tracker, the transverse and 
longitudinal impact parameters of the tracker track with respect to the event primary vertex, 
and the isolation variable. The muon isolation variable ( R ^ so ) is defined in a similar manner to 
that for electrons, but with a cone of radius A R = 0.4. Isolated muons must satisfy R^ so < 0.12. 
After all these requirements, muons are selected if they satisfy pj > 30 GeV and \rj\ < 2.1. 

Both electron and muon identification efficiencies are measured via the tag-and-probe tech¬ 
nique using inclusive samples of Z —> £ + £ events from data and simulation. Correction fac¬ 
tors are used to account for the difference in performance of the lepton identification between 
data and simulation. 
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5 Event selection 


Jets are reconstructed from PF candidates that are clustered with the anti-fcj algorithm It33l 
with a distance parameter of 0.5, using the FASTjET package |34| . Jet energy scale corrections 
obtained from data and simulation are applied to account for the response function of the 
combined calorimetry to hadronic showers and pileup effects 113511361. The jet pj resolution in 
simulation is adjusted to match that measured in data 1371 . Jet candidates are required to have 
pj > 30 GeV and ;/| < 4.0, and to satisfy a very loose set of quality criteria 11371 . The combined 
secondary vertex (CSV) b-tagging algorithm Il38l is used to identify jets from the hadronization 
of b quarks. The CSV algorithm exploits the large impact parameters and probable presence of 
a displaced vertex which are common in b-quark-initiated jets. This information is combined 
in a likelihood discriminant providing a continuous output between 0 and 1. In this search, a 
selected jet is considered to be b-tagged if it has a CSV discriminant value greater than 0.679 
and 1 p\ < 2.4. The b-tagging efficiency is approximately 70% (20%) for jets originating from 
a b (c) quark and the mistagging probability for jets originating from light quarks or gluons 
is approximately 2%. An event-by-event correction factor is applied to simulated events to 
account for the difference in performance of the b-tagging between data and simulation 11391 . 

Missing transverse energy (E™ ISS ) is measured as the magnitude of the vectorial pj sum of all 
PF candidates, taking into account the jet energy corrections. 


5 Event selection 

In semileptonic tt decays, two b quarks and two light quarks are produced. Therefore most of 
the selected signal events contain at least four jets. However, we set the requirement to be three 
or more rather than four or more identified jets in an event, since this is found to improve the 
search sensitivity by 10%. In addition, we require at least one b-tagged jet ("b jet") in the event, 
and only one identified isolated lepton. 

Signal events usually have larger Ey ,ss than the backgrounds because of two DM particles, 
neither of which leave any energy in the detector. Events are therefore required to have Ey ,ss > 
160 GeV. These selection criteria are referred to as the "preselection". After preselection, the 
dominant backgrounds are from tt and W+jets production. Other backgrounds include single 
top, Drell-Yan and diboson production. The QCD multijet contribution to the background is 
negligible because of the requirements of a high-p T isolated lepton, large E™ 1SS , and a b-tagged 
jet. 

To improve the search sensitivity, we further select events with Ey ISS > 320 GeV. The remaining 
W+jets and most tt backgrounds contain a single leptonically decaying W boson. The trans¬ 
verse mass, defined as Mj = \/ 2E™ iss py (1 — cos(A (p)), where p' T is the transverse momentum 
of the lepton and A (p is the opening angle in azimuth between the lepton and p| 11iss vector, is 
constrained kinematically to Mj < Mw for the on-shell W boson decay in the tt and W+jets 
events. For signal events, off-shell W boson decays, and tt dilepton decay channel, Mj can ex¬ 
ceed Mw- Therefore a requirement of Mj > 160 GeV is applied to increase the discrimination 
of the background relative to the signal. 

The dominant background with large Mj arises from dileptonic tt events where one of the lep¬ 
tons is unobserved, illustrated in Fig. [ 2 ] The variable Il40l is exploited to further reduce this 
type of background. This variable is defined as the minimal "parent" particle mass compati¬ 
ble with all the transverse momenta and mass-shell constraints, assuming two identical parent 
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particles, each of mass m y , decaying to bW: 


= min 


nu 


consistent with: 


P 1 


Pi = pf ss , rl = o, (pi + Pt ) 2 = vl = m: 


W' 


(pi + P€ + Pbl) 2 = (P2 + Pb2) 2 = m 2 


(3) 


where the momentum of the W boson that decays to an unreconstructed lepton is indicated 
by p 2 , and the momentum of the neutrino from the decay of the other W boson is indicated 
by pi. In particular, the intermediate W bosons are assumed to be on-shell, thus adding more 
kinematic information to suppress dileptonic tt events where one lepton is lost. In tt events, the 
distribution has a kinematic end-point at the top-quark mass, assuming perfect measure¬ 
ments with the detector. By contrast, this is not the case for signal events where two additional 
DM particles are present. The calculation of requires that at least two b jets be identified 
and be paired correctly to the lepton. When only one b jet is selected, each of the first three re¬ 
maining highest px jets is considered as the second b jet. When two or more b jets are selected, 
all the b jets in the event are used. The value is then calculated for all possible jet-lepton 
combinations and the minimum value is taken as the event discriminant. We select events with 
M* > 200 GeV. 



Figure 2: Schematic of a dileptonic tt event where only one lepton is reconstructed ll4C)ll . This 
represents the dominant type of tt background to this search. The momentum of the W boson 
that decays to an unreconstructed lepton is indicated by p 2 , and the momentum of the neutrino 
from the decay of the other W boson is indicated by pi. The same notation is used in Eq. (|3j) . 


In addition, the jets and the p™ iss tend to be more separated in cp in signal events than in tt 
background. We therefore require the minimum opening angle in (p between each of the first 
two leading jets and pj liss to be larger than 1.2. In summary, the signal region (SR) for our 
search is E™' ss > 320 GeV, Mj > 160 GeV, > 200 GeV and min Ac /?(71 2 , pj' ss ) > 1.2. These 
selection criteria are optimized based on the expected significance for DM masses between 1 
and 1000 GeV. 

Figure |3] shows the distributions of E™ ISS , Mj, M^, and min Acpl/i^, p™' ss ) after applying all 
other selections except the one plotted, indicating their power of discrimination between signal 
and background. In these distributions, the tt+jets and W+jets backgrounds have been adjusted 
by the scale factors (SF), as described in section [ 6 ] 
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5 Event selection 


19.7 fb 1 (8 TeV) 


19.7fb' 1 (8TeV) 
Drell-Yan 
Di boson 
M z = 1 GeV 
M, = 600 GeV 



300 
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Figure 3: Distributions of Ey ISS , My, M^, and min A(p(j^ 2 , p™ iss ) after applying SFs for tt+jets 
and W+jets backgrounds, as described in section [6] Each distribution is plotted after apply¬ 
ing all other selections, which are indicated by the arrows on the relevant distributions. Two 
simulated DM signals with mass M x of 1 and 600 GeV and an interaction scale M* of 100 GeV 
are included for comparison. The hatched region represents the total uncertainty in the back¬ 
ground prediction. The last bin of the E™ ISS , My and distributions includes the overflow. 
The horizontal bar on each data point indicates the width of the bin. 
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6 Background estimation 

Standard model backgrounds are estimated from simulation, with data-to-simulation SFs ap¬ 
plied to the dominant backgrounds from tt+jets and W+jets. 

Two control regions (CR) are defined to extract these SFs. One is the preselection with the 
additional requirement of Mj > 160 GeV (CR1). The sample in CR1 is dominated by tt+jets 
background. The other (CR2) is defined the same way as CR1 except that no jet should satisfies 
the b-tag requirement, resulting in a sample enriched in W+jets events. The subdominant back¬ 
grounds are subtracted from the distributions observed in data in order to obtain a data sample 
that has only tt+jets and W+jets background contributions. The tt+jets and W+jets SFs are then 
obtained by matching simultaneously to data the Mj distribution in CR1 and the E™ ISS distribu¬ 
tion in CR2. The obtained SFs for tt+jets and W+jets are 1.11 ± 0.02 (stat) and 1.26 ± 0.06 (stat), 
respectively. These SFs are propagated to the SR to estimate the background. The level of DM 
signal contamination in the two CRs is estimated to be small and therefore has negligible im¬ 
pact on the background estimation in the SR. Figures [ 4 ] and [ 5 ] show the distributions of Ej' ss , 
Mj, Mj 2 , and min 2 , p™ iss ) with the SFs applied in CR1 and CR2, respectively. The data 
are in good agreement with expectations from SM background. 


7 Systematic uncertainties 

The normalization and shape of the distributions used to establish a possible DM signal are 
subject both to experimental and theoretical uncertainties. 

The data-to-simulation SFs for tt+jets and W+jets are extracted from the CRs, as described 
in the previous section. For the background estimation, the use of SFs largely removes the 
uncertainties from the integrated luminosity, lepton identification and trigger efficiencies, and 
from cross sections of the two backgrounds. Other systematic uncertainties can be constrained 
by refitting the data in the CRs, as described in the following. 

The tt+jets and W+jets SFs are obtained from CRs in which other backgrounds are present as 
well. We conservatively assign a 50% uncertainty for other backgrounds to account for possible 
missing higher order terms as well as mismodelling of kinematic properties from the simula¬ 
tion. This uncertainty results in a change of 5% and 9% for the tt+jets and W+jets SFs, respec¬ 
tively. Propagating these changes to the SR, the impact on the total background prediction is 
found to be 10%. 

The stability of the SFs is checked through changes in the definitions of the CRs. These include 
tightening the E™ 1SS requirement or applying selections on Mjj, and min Ac/)f/i 2 , p™' ss )- An 
uncertainty of 40% for the W+jets SF is assigned from these CR tests. No significant change is 
observed in the SF for tt+jets. 

The Pj distributions of top quarks in the tt+jets simulation is reweighted to match the data. The 
reweighting uncertainty is estimated by changing the nominal reweighting factor to unity or to 
the square of the reweighting factor, resulting in a change of ±14% for the tt+jets SF and only 
negligible impact on the W+jet SF. Propagating these SFs to the SR, a systematic uncertainty of 
10% is estimated for the tt+jets background prediction from the reweighting. The stability of 
the tt+jets background prediction is also checked by varying the MadGraph factorization and 
renormalization scale parameters, or the scale parameter for the matrix element and parton 
shower matching, by a factor of two. The resulting predictions are consistent with the nominal 
tt+jets background prediction. 
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7 Systematic uncertainties 
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Figure 4: Distributions of Ej' ss , Mj, and min A(p(j\ r 2 , pT" ss ) in CR1 after applying the SFs 
for tt+jets and W+jets backgrounds, as described in section pi Two simulated DM signals with 
mass M x of 1 and 600 GeV and an interaction scale AT, of 100 GeV are included for comparison. 
The hatched region represents the total uncertainty in the background prediction. The error 
bars on the data-to-background ratio take into account both the statistical uncertainty in data 
and the total uncertainty in the background prediction. The last bin of the E™ ISS , Mj, and Mj 2 
distributions includes the overflow. The horizontal bar on each data point indicates the width 
of the bin. 
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Figure 5: Distributions of Ej' ss , Mj, and min At/T/i^, pT" ss ) in CR2 after applying the SFs 
for tt+jets and W+jets backgrounds, as described in section pi Two simulated DM signals with 
mass M x of 1 and 600 GeV and an interaction scale AT, of 100 GeV are included for comparison. 
The hatched region represents the total uncertainty in the background prediction. The error 
bars on the data-to-background ratio take into account both the statistical uncertainty in data 
and the total uncertainty in the background prediction. The last bin of the E™ ISS , Mj, and Mj 2 
distributions includes the overflow. The horizontal bar on each data point indicates the width 
of the bin. 
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8 Results 


The remaining dominant experimental systematic uncertainties are from corrections in jet en¬ 
ergy scale and resolution. Correction factors are separately varied by ±1 standard deviation 
and Ej' ss is recalculated accordingly These changes in the jet energy scale and resolution 
correction factors contribute uncertainties of 4% and 3% in the estimate of the background, re¬ 
spectively The uncertainties in the background yield due to b-tagging correction factors are 
estimated to be 1.0% and 1.8% for heavy-flavour and light-flavour jets, respectively. The uncer¬ 
tainty in the pileup model contributes an uncertainty of 2.0% in the background estimate. 

The theoretical uncertainty related to the choice of the PDF set is evaluated by reweighting the 
background samples using three PDF sets: CT10 ED . MWST2008 |[42l , and NNPDF2.3 11431 . 
following the PDF4LHC recommendation Il44l 1451. For each PDF set, an uncertainty band is 
derived from the different error PDF sets, including the uncertainties due to the strong coupling 
constant ag. The envelope of these three error bands is taken as the PDF uncertainty, which 
leads to a 2.6% uncertainty in the background estimate. 

Table[l]summarizes the systematic uncertainties and their impact on the background prediction 
in the SR. 

The following sources of systematic uncertainty associated with the signal expectation are 
taken into account. The integrated luminosity is measured with precision of 2.6% Ii46l . Lepton 
trigger and identification efficiencies are measured with a precision of 2% and 1%, respectively. 
Uncertainties in the jet energy scale and resolution correction factors yield uncertainties of 2- 
3% and less than 1%, respectively, depending on the mass hypotheses for the DM particle. 
Uncertainties in the b-tagging correction factors for heavy-flavour and light-flavour jets yield 
uncertainties of 3-4% and less than 1%, respectively. 

Table 1: Systematic uncertainties from various sources and their impact on the total background 
prediction. 


Source of systematic uncertainties 

Relative uncertainty on 
total background (%) 

50% normalization uncert. of other bkg in deriving SFs 

10 

SFw+jets (CR tests) 

13 

tt+jets top-quark pj reweighting 

3.9 

Jet energy scale 

4.0 

Jet energy resolution 

3.0 

b-tagging correction factor (heavy flavour) 

1.0 

b-tagging correction factor (light flavour) 

1.8 

Pileup model 

2.0 

PDF 

2.6 


8 Results 

Table [2] lists the number of events observed in the SR, along with the background prediction 
and expected number of signal events for a DM particle with mass of M x = 1 GeV and an in¬ 
teraction scale M* = 100 GeV. We observe no excess of events in the SR and set 90% confidence 
level (CL) upper limits on the production cross section of DM particles in association with a 
pair of top quarks. The choice of 90% CL is made in order to allow direct comparisons with 
related limits from astrophysical observations. A modified-frequentist CL S method |47, 48J is 
used to evaluate the upper limits, with both statistical and systematic uncertainties taken into 
account in the limit setting. 
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Table 2: Expected number of background events in the SR, expected number of signal events 
for a DM particle with the mass M x = 1 GeV, assuming an interaction scale AT, = 100 GeV, and 
observed data. The statistical and systematic uncertainties are given on the expected yields. 


Source 

Yield (±stat ±syst) 

tt 

8.2 ±0.6 ±1.9 

W 

5.2 ±1.8 ±2.1 

Single top 

2.3 ±1.1 ±1.1 

Diboson 

0.5 ±0.2 ±0.2 

Drell-Yan 

0.3 ±0.3 ±0.1 

Total Bkg 

16.4 ±2.2 ±2.9 

Data 

18 


Table [3] shows the signal efficiencies and the observed and expected upper limits on the pp —> 
tt + XX production cross section for seven mass hypotheses of the DM particle. The relatively 
low values of signal efficiencies of 1-3% are mostly due to the requirement of E™ ISS > 320 GeV. 
Cross sections larger than 20 to 55 fb are excluded at 90% CL for DM particles with mass rang¬ 
ing from 1 to 1000 GeV. Interpreting the results in the context of a scalar interaction between 
DM particles and top quarks, we set lower limits on the interaction scale M*, shown in Fig. [6] 
Assuming a DM particle with a mass of 100 GeV, values of the interaction scale below 118 GeV 
are excluded at 90% CL. 


Table 3: Expected number of signal events in SR assuming an interaction scale AT = 100 GeV, 
signal efficiencies, and observed and expected limits at 90% CL on production cross sections 
for pp —» tt + XX/ f° r various DM particle masses. 


M x (GeV) 

Yield (±stat ±syst) 

Signal efficiency (%) (±stat ±syst) 

^(fb) 

( fc ) 

1 

38.3 ±0.7 ±2.1 

1.01 ±0.02 ±0.05 

47+ 21 
^ —13 

55 

10 

37.8 ±0.7 ±2.1 

1.01 ±0.02 ±0.05 

46+ 21 

54 

50 

35.1 ±0.6 ±1.9 

1.20 ±0.02 ±0.06 

39+^ 

45 

100 

30.1 ±0.4 ±1.7 

1.46 ±0.02 ±0.07 

32± 44 

37 

200 

18.0 ±0.2 ±1.0 

1.73 ±0.02 ±0.08 

27+1 2 

32 

600 

1.26 ±0.02 ±0.07 

2.40 ±0.03 ±0.11 

19 +9 

23 

1000 

0.062 ± 0.001 ± 0.003 

2.76 ±0.04 ±0.13 

17 +8 

20 


As shown in Eq. (JTJ), DM production is modeled by an ELT, an approximation that has some 
important limitations. Firstly, the EFT approximation is only valid when the momentum trans¬ 
fer Q tr is small compared to the mediator mass. Secondly, the couplings should not exceed 
the perturbative limit. Unfortunately, both of these conditions depend on the details of the 
unknown new physics being approximated by the EFT. For example, if we consider a model 
with s-channel exchange between the top quarks and the DM particles and a coupling equal 
to the perturbative limit g = y/g x gt = 47r, where g x and are the coupling constants of the 
mediator to DM particles and top quarks, respectively, then we can derive a lower bound on 
M*, \JMl/m t > M x /2tz, where m t is the mass of the top quark [|3l!49|. The region of parameter 
space in the exclusion plane that does not meet the perturbative condition for the validity of 
the EFT is indicated by the hatched area in Fig. [6] 

In addition to this minimal requirement, we also test the validity of the EFT approximation 
with respect to the momentum transfer condition. For the same s-channel mediator scenario, 
Qtr is estimated as the invariant mass of two DM particles ( M xx ) as shown in Fig. [7] The 
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9 Summary 



Figure 6: Observed exclusion limits in the plane of DM particle mass and interaction scale, with 
the region below the solid curve excluded at a 90% CL. The background-only expectations are 
represented by their median (dashed line) and by the 68% and 95% CL bands. A lower bound 
of the validity of the EFT is indicated by the upper edge of the hatched area. The four curves, 
corresponding to different g and R values, represent the lower bound on M, for which 50% and 
80% of signal events have a pair of DM particles with an invariant mass less than g\J M*/m t , 
where g = 4n and g = 2n respectively. These curves indicate further restrictions on the 
applicability of EFT, as explained in the text. 


EFT approximation is then valid if M xx < g\JM\Jni x The fraction of simulated signal events 
that satisfy this requirement (R) is reported for given values of g and M*. For g = 4 k and 
g = 2 tc, contours are overlaid in Fig. [6]that indicate where in the exclusion plane 50% or 80% of 
simulated signal events passing the analysis selection criteria satisfy the momentum transfer 
condition. If instead of drawing such a contour we fix M* at the 90% CL lower limit obtained in 
this analysis, then 89% (46%) of simulated signal events passing the analysis selection criteria 
satisfy the momentum requirement for g = 4n(2n) and M x = lGeV. These fractions drop 
to 63% (5%) for M x = 200 GeV. No simulated signal events passing the analysis selection 
criteria are found to satisfy this requirement for M x > 600 GeV. For these reasons, the 90% 
CL constraints on A/L obtained in this analysis cannot be considered generally applicable, but 
should only be interpreted in models with large DM coupling. 

The limits on the interaction scale M* can be translated to limits on the DM-nucleon scatter¬ 
ing cross section |3). Figure [8] shows the observed 90% CL upper limits on the DM-nucleon 
cross section as a function of the DM mass for the scalar operator considered in this paper. 
More stringent limits are obtained relative to current direct DM searches in the mass region of 
less than ~6 GeV. In this region, DM-nucleon cross sections larger than 1-2 x 10 42 cm 2 are 
excluded. 


9 Summary 

A search has been presented for the production of dark matter particles in association with top 
quarks in single-lepton events with the CMS detector at the LHC, using proton-proton collision 
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19.7 fb 1 (8 TeV) 



Invariant mass of DM pairs M (GeV) 

XX 

Figure 7: Invariant mass of two DM particles M AA in selected signal events, for several DM 
mass hypotheses. 

data recorded at y/s = 8 TeV and corresponding to an integrated luminosity of 19.7 fb -1 . No ex¬ 
cess of events above the SM expectation is found and cross section upper limits on this process 
are set. Cross sections larger than 20 to 55 fb are excluded at 90% CL for dark matter particles 
with the masses ranging from 1 to 1000 GeV. Interpreting the findings in the context of a scalar 
interaction between dark matter particles and top quarks in the framework of an effective field 
theory, lower limits on the interaction scale are set. Assuming a dark matter particle with a 
mass of 100 GeV, values of the interaction scale below 118 GeV are excluded at 90% CL. These 
limits on the interaction scale are comparable to those obtained from a similar search by the 
ATLAS Collaboration Ifl2l . In the case of an s-channel mediator, they are only valid for large 
values of the coupling constant, where the effective field theory approximation holds for most 
signal events. These limits are interpreted as limits on the dark matter-nucleon scattering cross 
sections for the spin-independent scalar operator. For dark matter particles with masses below 
6 GeV, more stringent limits are obtained from this search than from direct dark matter detec¬ 
tion searches. Dark matter-nucleon cross sections larger than 1-2 x 10 42 cm 2 are excluded at 
90% CL. 
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19.7 fb 1 (8 TeV) 



Figure 8: The 90% CL upper limits on the DM-nucleon spin-independent scattering cross sec¬ 
tion (ct^ 1 n ) as a function of the DM particle mass for the scalar operator considered in this paper. 
Also shown are 90% CL limits from various direct DM search experiments 15014541 . 
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